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I '■ Abstract 
(N ; 

. We present results on the determination of the observable ratio R = BR{H^ — > 

\ T^u^) / BR{H^ tb) of charged Higgs boson decay rates as a discriminant quantity 

' between Supersymmetric and non-Supersymmetric models. Simulation of measure- 

r~| ' ments of this quantity through the analysis of the charged Higgs production process 

Q-i! gb tbH^ and relative backgrounds in the two above decay channels has been per- 

Q^' formed in the context of ATLAS. A ~ 12 — 14% accuracy on R can be achieved for 

^ ■ tan/3 = 50,mfj± = 300 — 500 GeV and after an integrated luminosity of 300 fb~^. 

• • , With this precision measurement, the Large Hadron Collider (LHC) can easily dis- 

' criminate between models for the two above scenarios, so long as tan/3 > 20. 

' In this note we investigate the production of charged Higgs bosons in association with 
top quarks at the LHC, from the experimental and theoretical point of view, by studying 
hadronic {H~^ — >■ tb) and leptonic {H^ r+z/) decay signatures. The interest of this 
investigation is many-fold. 

• The discovery of a charged Higgs boson will point immediately to the existence of 
some extension of the Standard Model (SM). 

• The associated production of a charged Higgs boson with a top quark {pp — > H~^t + 
X) [1,2] is only relevant at large values of tan/?^, a regime where Higgs boson ob- 
servables receive large Supersymmetric (SUSY) radiative corrections. 



*Electronic addresses: Ketevi.AdikleAssamagan@cern.ch, jaume.guasch@psi.ch, stefano@hep.phys. 
soton. ac . uk , siannah@mppmu . mpg . de . 

^It is in principle also relevant at very low values of tan/3 (say, < 1). In practise, this tan/3 regime is 
excluded in the Minimal Supersymmetric Standard Model (MSSM) from the negative neutral Higgs search 
at LEP [3]. Hence, hereafter, we will refrain from investigating the low tan/3 case. 
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• While SUSY radiative effects might be difficult to discern in the production cross- 
sections separately, they will appear neatly in the following relation between the two 
mentioned channels: 

^ _ ajpp H+t + X r+z/t + X) 

a{pp H+i +X ^tbi+X) ' ^ ^ 

• In fact, in the ratio of (P), the dependence on the production mode (and on its 
large sources of uncertainty deriving from parton luminosity, unknown QCD radiative 
corrections, scale choices, etc.) cancels out: 

^ ^ BRiH+ T+u) ^ TiH+ ^ r+z.) 

BR{H+ tb) r{H+ -^tb) ' ^ ' 

From these remarks it is clear that the quantity R is extremely interesting both experi- 
mentally and theoretically in investigating the nature of Electro- Weak Symmetry Breaking 
(EWSB). 

In the MSSM, Higgs boson couplings to down-type fermions receive large quantum 
corrections, enhanced by tan/5. These corrections have been resummed to all orders in 
perturbation theory with the help of the effective Lagrangian formalism for the tbH^ 
vertex [4,5]. The 6-quark Yukawa coupling, /z^ , is related to the corresponding running 
mass at tree level by /i;, = mi,/vi. Once radiative corrections are taken into account, due 
to the breaking of SUSY, this relation is modified to mb = h^vi (1 -|- Am;,) [4,5], where 
Vi is the Vacuum Expectation Value (VEV) of the Higgs doublet Hi and Am^, is a non- 
decoupling quantity that encodes the leading higher order effects. Similarly to the 6-quark 
case, the relation between m,- and the r-lepton Yukawa coupling, h.^, is also modified by 
quantum corrections. Am,-. We adopt in our analysis the effective Lagrangian approach 
by relating the fermion mass to the Yukawa coupling via a generic Arrij- (/ = 6, r), 

hf = ^-^ = ^-^ {v = {vl + vlfl\ tan/? = ^), (3) 
^ vi 1 + Am^ fcos/?l + Am^ ^ ' V2 

in which the resummation of all possible tan/? enhanced corrections of the type (a^tan/?)"' 
is included [4,5]. The leading part of the (potentially) non-decoupling contributions pro- 
portional to soft-SUSY-breaking trilinear scalar couplings {Af) can be absorbed in the 
definition of the effective Yukawa coupling at low energies and only subleading effects sur- 
vive [5]. Therefore, the expression (jHl) contains all (potentially) large leading radiative 
effects. The SUSY-QCD contributions to Am^ are proportional to the Higgsino mass pa- 
rameter Arrtf, ~ /i, while the leading SUSY-EW contributions behave like Am^, ~ ^At [6]. 
Thus, they can either enhance or screen each other, depending on the sign of At. It is 
precisely these effects that will allow us to distinguish between different Higgs mechanisms 
of EWSB. For example, the analysis of these corrections in the ratio of neutral Higgs boson 
decay rates, R' = BR{H bb)/BR{H t~^t~), revealed large deviations from the SM 
values for several MSSM parameter combinations [6]. Extensive theoretical analyses of 
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one-loop corrections to both neutral and charged Higgs boson decays have been performed 
in [6-17]. We now explore the one- loop MSSM contributions to the ratio of the branching 
ratios (BRs) of a charged Higgs boson in (j2|), which at leading order (and neglecting 
kinematical factors) is given hjR = h'^/Shl in the large tan/3 limit. The SUSY corrections 
to the H~^tb vertex entering the decay processes t — *■ H~^b and th have been anal- 

ysed in [7-10], where it was shown that they change significantly the Tevatron limits on 
mH± [10]. They were further explored in the production process pp{pp) H^tb at LHC 
and Tevatron in [11-13,18], where they were shown to shift significantly the prospects for 
discovery of a charged Higgs boson at both colliders. 

Here, we have performed a detailed phenomenological analysis for the LHC of charged 
Higgs boson signatures, by using the subprocess gb — >• H^i. The QCD corrections to 
this channel are known to next-to-leading (NLO) [19]. However, we have normalised our 
production cross-section to the LO result, for consistency with the tree-level treatment of 
the backgrounds^. In our simulation, we have let the top quarks decay through the SM-like 
channel t — > W~^b. In the hadronic decay channel of the charged Higgs boson {H~^ —>■ tb) we 
require one of the two W^s emerging from the decay chain H^t — > {tb)t — > {W^b)b{bW~) 
to decay leptonically, to provide an efficient trigger, while the other W is forced to decay 
hadronically, since this mode provides the largest rate and in order to avoid excessive 
missing energy. The r-lepton in the r^u decay mode is searched for through 

hadronic one- and multi-prong channels. In summary, the experimental signatures of the 
two production channels under investigation are {I = e,fi): 

pp{gb) ^ Hn^{r+u)i^r+u{jfb) , (4) 
pp{gb) ^ Hn^{tb)i^{jj[lu]b)b{lu[jj]b) . (5) 

(In the numerical analysis we always combine the signals in (j3]) and ^ with their charged- 
conjugated modes.) 

The Monte Carlo (MC) simulation has been performed using PYTHIA (v6.217) [20] for 
the signal and most of the background processes. (We have cross-checked the signal cross- 
section with [19].) We have used HDECAY [21] for the Higgs boson decay rates. One of the 
background processes (the single-top one: see below) has been generated with TopRex [22] 
with a custom interface to PYTHIA. We have used ATLFAST [23] for the detector simulation. 
(Further details of the detector can be found in [24,25].) We have adopted the CTEQ5L [26] 
parton distribution functions in their default PYTHIA v6.217 setup and we have used running 
quark masses derived from the pole values mf°'*' = 175 GeV and ml"^'^ = 4.62 GeV. The 
TAUOLA [27-29] package was interfaced to the PYTHIA event generator for treatment of the 
r-lepton polarisation. 

The leptonic decay channel of the charged Higgs boson provides the best probe for the 
detection of such a state at the LHC. In fact, it turns out that despite the small branching 

^In all the analysis we disregard the subleading QCD and SUSY corrections which affect the signal and 
the background, and will take only into account the leading SUSY corrections to the signal cross-section, 
which are absent in the background processes. 
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mH± = 350 


mH± = 500 


tt 


WH 


ax BR 


99.9 fb 


30.7fb 


79.1 pb 


16.3 pb 


Events 


29958 


9219 


2.3 X 10*^ 


4.89 X 10^ 


Events after cuts 


174 


96 


17 


3 


Efficiency 


0.6% 


1% 


8 X 10-"^ 


6 X IQ-^ 


S/B 


7.9 


4.4 




s/Vb 


37.1 


20.5 


Poisson 


23.1 


14.6 



Table 1: Tlie signal and background cross-sections, the number of events before cuts, the 
number of events after all cuts, the total efficiency, the signal-to-background ratios (S/B), 
and the signal significances (Gaussian and Poisson) for the detection of the charged Higgs 
in the riy channel at the LHC, for 300 fb~^ integrated luminosity and tan/3 = 50. 

ratio BR[H~^ t^^^), the r-lepton affords an efficient trigger to observe this channel. The 
production rates a x BR{H^ r+z/) x BR{W jj) are shown in Tab. ^ The main 
background processes in this channel are: top-pair production with one of the W's decaying 
into Tu {gg tt ^ jj brub) and W^t associated production {gb W^i T^vt). 

We have used the following trigger conditions: hadronic r-jet (p^- > 30 GeV); a 6-tagged 
jet (py > 30 GeV) and at least two light jets (p^ > 30 GeV). We apply afterwards a 6-jet 
veto to reject the ti QCD background. As there is no isolated lepton (electron or muon) in 
the final state, the observation of this channel requires a multi-jet trigger with a r-trigger. 
After reconstructing the jet-jet invariant mass mjj and retaining the candidates consistent 
with the W-hosow mass, \mw — "^jjl < 25 GeV, the jet four-momenta are rescaled and the 
associated top quark is reconstructed by minimising = i^jjb ~ ^t)"^- Subsequently, a 
sufficiently high threshold on the pt of the r-jet is required, > 100 GeV. The background 
events satisfying this cut need a large boost from the H^-boson. This results in a small 
azimuthal opening angle A0 between the r-jet and the missing transverse momentum, 
For background suppression we then have the cut A(j){'^j,,p'!p) > 1. Besides, the missing 
transverse momentum is harder for the signal than for the background while the differences 
between their distributions in azimuthal angle and missing transverse momentum increase 
with increasing mH±. These effects are well cumulated in the transverse mass, ttlt = 
>y2p'Jp'^j, [1 — cos(A0)], which provides good discrimination between the signal and the 
backgrounds, as shown in Fig.Q (Further details of this kind of studies are available in [24, 
25].) The discussed set of cuts reduces the total background by six orders of magnitude 
while the signal is only suppressed by two orders. The production rates and total detection 
efficiency (including detector acceptance, b- and r-identification, pileup and the effect of 
cuts) are also shown in Tab. ^for an integrated luminosity of 300 fb~^. We can see that 
the signal rates are large enough to indeed consider tu a. golden channel for the 

discovery at large tan/3. 

The production rates a x BR{H^ tb) x BRCW^W^ — » jjl^) are shown in Tab. El 
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Figure 1: Transverse mass itlt distribution for signal and total background taking into 
account the polarisation of the r-lepton, for an integrated luminosity of 30fb~^. A final 
cut TTiT > 200 GeV was used for the calculations of the signal-to-background ratios and 
for the signal significances. 

The decay mode —>■ tb has large QCD backgrounds at hadron colliders that come 
from ttq production with ti — > WbWb —>■ Ivbjjb. However, the possibility of efficient b- 
tagging has considerably improved the situation [30-32] . We search for an isolated lepton 
(pj. > 20 GeV, > 8 GeV), three 6-tagged jets (p^ > 30 GeV) and at least two non 6-jets 
{pip > 30 GeV). We retain the jet-jet combinations whose invariant masses are consistent 
with the VT-boson mass, \mw — mjj\ < 25 GeV, then we use the W^-boson mass constraint 
to find the longitudinal component of the neutrino momentum in lu, by assuming 

that the missing transverse momentum belongs only to the neutrino. Subsequently, the two 
top quarks entering the H~^i (tb)i {jj[lh']b)b{lh'[jj]b) decay chain are reconstructed, 
retaining the pairing whose invariant masses mi^b and mjjb minimise = {^t — fnivbY + 
{mt — rrijjh)^. The remaining 6-jet can be paired with either top quark to give two charged 
Higgs candidates, one of which leads to a combinatorial background. The expected rates 
for signal and background (after the mentioned decays) are shown in Tab.|21 (This analysis 
is presented extensively in [24,25].) 

At this point, we have two charged Higgs candidates: tifes or t2&3- Assuming that 
the charged Higgs is discovered through the rv channel and its mass determined 

from the tv transverse mass distribution [24, 25] , the correct charged Higgs candidate in 
the tb channel can be selected by using the measured m/^± as a constraint. This is done 
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gb H^t Ivbbjjb 
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Figure 2: The signal and the background distributions for the reconstructed invariant mass 
rritb of = 350 GeV, tan /3 = 50 and an integrated luminosity of 30 fb~^. Assuming that 
the charged Higgs is discovered in the tu channel, one can use as a constraint 

to reduce the combinatorial background: this is shown on the right plots. 



by selecting the candidate whose invariant mass is closest to the measured charged Higgs 
mass: = (^tfe ~ ^h±Y- The signal distribution for the reconstructed invariant mass 
TJitb for a charged Higgs boson weighing 350 GeV, with tan /3 = 50 and after integrated 
luminosity of 30fb~^, is shown in Fig. |21 We can see that for the tb decay some 

irreducible combinatorial noise still appears even when the mj^± constraint is included. 
In addition, for the background, we have found that the mj^± constraint reshapes the 
distributions in gg — > ttX in such a way that no improvement in the signal-to-background 
ratio and signal significance is further observed. Finally, recall that the knowledge of the 
shape and the normalisation of the reshaped background would be necessary for the signal 
extraction. For these reasons, we did not use the mH± constraint for the results shown in 
this work. The subtraction of the background can then be done by fitting the side bands 
and extrapolating in the signal region which will be known from the mH± determination 
in the —>■ tv channel: however, this would be possible only for Higgs masses above 
300 GeV - see Fig. |21 The signal and background results are summarised in Tab. |21at an 
integrated luminosity of 300 fb~^ for different values of mH± and tan/5 = 50. It is shown 
that it is difficult to observe signals in this channel above ~ 400 GeV, even with the 
mH± constraint. For masses above ^ 400 GeV the signal significance can be enhanced 
by using the kinematics of the three-body production process gg — > H'^ib [11-13,30-32]. 
We assume a theoretical uncertainty of 5% on the branching ratios, BR s. Previous AT- 



6 





mu± = 350 


mH± = 500 


ttq 

J. 


ax BR 


248.4 fb 


88 fb 


85 pb 


Events 


74510 


26389 


2.55 X 10^ 


Events after cuts 


2100 


784 


59688 


Efficiency 


2.8% 


3% 


0.2% 


S/B 


0.035 


0.013 




s/Vb 


8.6 


3.2 



Table 2: The signal and background cross-sections, the number of events before cuts, 
the number of events after all cuts, total efficiency, S/B, and the signal significances for 
the detection of the charged Higgs in the th channel at the LHC, for 300 fb~^ integrated 
luminosity and tan/? = 50. 

LAS studies have shown the residual gg tt shape and normalisation can be determined 
to 5% [24,25]. The scale uncertainties on jet and lepton energies are expected to be of 
the order 1% and 0.1% respectively [24,25]. As explained above, for mH± > 300 GeV, the 
side band procedure can be used the subtract the residual background under the — >■ th 
signal: we assume also a 5% uncertainty in the background subtraction method. Thus, 
the statistical uncertainties can be estimated as \Jl/ S. The uncertainty in the ratio R are 
dominated by the reduced knowledge of the background shape and rate in the —>■ th 
channel. The cumulative results for the two channels are summarised in Tab. [Slat an inte- 
grated luminosity of 300 fb~^. Here, the final result for the ratio R is obtained by correcting 
the visible production rates after cuts for the total detection efficiency in Tabs. [T] and El and 
by the decay BR s of the W^-bosons. The simulation shows that the above mentioned ratio 
can be measured with an accuracy of ~ 12 — 14% for tan (3 = 50, for mH± = 300 — 500 GeV 
and at an integrated luminosity of 300 fb~^. 

We turn now to the impact of the SUSY radiative corrections. Their role is twofold. 
Firstly, by changing the value of the Yukawa coupling they change the value of the ob- 
servable R. Secondly, they can change the value of the production cross-section a{pp — > 
H'^t + X), hence shifting significantly (by as much as 100 GeV) the range of charged Higgs 
masses accessible at the LHC [11-13]. 

To explore the second consequence, we rely on the fact that the bulk of the SUSY 
corrections to the production cross-section is given by the Yukawa coupling redefinition 
in [11-13, 19]. By neglecting the kinematic effects, taking the large tan/5 limit and 
assuming that the dominant decay channel of the charged Higgs boson is — »• th (large 
mass limit), we can estimate the corrected production rates. For simplicity, we show only 
the contributions to Am^, since they are the dominant ones: 

a^^^igh^ H+i^tht) 

a^°"{gh H+t T+ut) 



a'°"igb ^ H+t) X BR'^'^iH^ ^ th) ~ "^"^^^ 7 ^"^P 



1 + Ami, 



7 





= 350 


mH± = 500 


Signals rv/th 


174 / 2100 = 0.08 


96 / 784 = 0.12 


Signals (corrected) rv/th 


0.18 


0.16 


Systematics unc. 


~ 9% 


~ 9% 


Total unc. 


12% 


14% 


Theory 


0.18 


0.16 



Table 3: Experimental determination of the ratio for 300 fb~^ and tan/5 = 50. Shown 
are: the signal after cuts, the signal after correcting for efficiencies and branching ratios, 
the systematic uncertainty, the total combined uncertainty, and the theoretical prediction 
(without SUSY corrections). 

a^igb H+i) T(H+ r+v) 

~ — — X = i' 

(l + AmJ2 V^H^ ^ th) ■ ^ 

This very simple exercise shows that the production rate in the r-channel is fairly inde- 
pendent of the SUSY radiative corrections and therefore the tree-level analysis performed 
above can (to a very good approximation) be used for our original purposes. Actually, 
once we take into account kinematical effects, the r-channel will receive small (negative) 
corrections in the low charged Higgs mass range. However, in this range, BR{H^ —>■ r^v) 
is quite large and one should not fear to loose the signal. Quite the opposite, the hadronic 
th production channel receives large radiative corrections. These corrections can be ei- 
ther positive (enhancing the signal, and therefore the significance in Tab. El) or negative 
(reducing it, possibly below observable levels)^. In Fig. O^i, we show the discussed en- 
hancement / suppression factors as a function of tan/? for mH± = 350 GeV and a SUSY 
mass spectrum defined as SPS4 of the Snowmass Points and Slopes in [33], but choos- 
ing different scenarios for the sign of fi and At"^. (It is worth noting that the production 
rate for the t6-channel can be enhanced by a factor larger than 3 in some SUSY scenar- 
ios, which would enhance significantly the corresponding signal thus overcoming the low 
signal-to-background ratio of this channel.) 

We now turn our attention to the observable under analysis. Fig. |2b shows the predic- 
tion for the ratio i? as a function of tan/5 for the SPS4 scenario with a charged Higgs mass 
of mH± = 350 GeV. The value of R only depends on m/^± through kinematical factors and 
the dependence is weak for mH± ~ 300 GeV. In this figure, we also show the experimental 
determination carried out as before and repeated for each SUSY setup. From Fig. it is 
clear that radiative SUSY effects are visible at the LHC at a large significance. In partic- 
ular, the fi < scenarios can easily be discriminated, while the /i > ones will be more 
difficult to establish, due to the lower signal rate of the hadronic channel. This feature then 
also allows for a measurement of the sign of the /j, parameter. In contrast, since radiative 
corrections are independent of the overall SUSY scale, the observable R cannot provide us 

■^Alternative analyses may permit the signal to be seen even in this unfavourable case [11-13]. 
^The SPS4 spectrum is affected by moderate SUSY radiative corrections. 
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Production rates SUSY factor 
SPS4 m„.=350 GeV 



H<0 A,<0 
|1>0 A,<0 
|1<0 A,>0 
H>0 A,>0 
' (7(TV)/(t''(TV) 

o(tb)/a°(tb) 



0.3 - 



no SUSY 

|1<0A,<0 

|1>0A,<0 

|1<0A,>0 

|1>0A,>0 


J 


















SPS4 m„+=350 GeV 
BR(H*->T*v)/BR(H*->tB) 





30 40 
tanP 



(b) 



Figure 3: a) Production rates enhancement /suppression factors for the r and the th chan- 
nels; b) the SUSY correction to the rate Q. Plots as functions of tan/? for mH± = 350 GeV 
and a SUSY spectrum as in SPS4, but for different choices of the signs of fi and Af. Shown 
is also the experimental determination for each scenario. 



with an estimation of the typical mass of SUSY particles. Nonetheless, the information 
obtained in other production channels (e.g., neutral Higgs bosons or SUSY particles direct 
production) can be used to perform precision tests of the MSSM. 

To summarise, we have used the observable R = '^(pp^-'^ J^"*"*^) = ^ 

' cr{pp^H+t+X^tbt+X) BR{H+ —*tb) 

to discriminate between SUSY and non-SUSY Higgs models. This quantity is a theoreti- 
cally clean observable. The experimental uncertainties that appear in this ratio have been 
analysed in details through detailed phenomenological simulations. In the MSSM, R is af- 
fected by quantum contributions that do not decouple even in the heavy SUSY mass limit. 
We have quantitatively shown that an LHC measurement of R can give clear evidence for 
or against the SUSY nature of charged Higgs bosons. 

Acknowledgements 

KAA, JG and SM wish to thank the Les Houches workshop organisation and participants 
for the warm atmosphere of the workshop. JG acknowledges financial support from the 
Les Houches workshop organisation and SM from The Royal Society (London, UK). 



References 

[1] J. F. Gunion, Phys. Lett. B322, 125-130 (1994), liep-ph/931220ll 

[2] V. D. Barger, R. J. N. Phillips and D. P. Roy, Phys. Lett. B324, 236-240 (1994), 
.hep-ph/9311372, 



9 



[3] The LEP Higgs Working Group, Searches for the neutral Higgs bosons of the MSSM: 
Prehminary combined results using LEP data collected at energies up to 209-GeV, 
|hep-ex /0107030( 

[4] M. Carena, D. Garcia, U. Nierste and C. E. M. Wagner, Nucl. Phys. B577, 88-120 
(2000),"hep-ph/9912516J 

[5] J. Guasch, P. Hafliger and M. Spira, Phys. Rev. D68, 115001 (2003), 
(hep-ph/030510l| 

[6] J. Guasch, W. HoUik and S. Peiiaranda, Phys. Lett. B515, 367-374 (2001), 
[hep-ph/0106027| 

[7] J. Guasch, R. A. Jimenez and J. Sola, Phys. Lett. B360, 47-56 (1995), 
hep-ph/9507461 . 

[8] J. A. Coarasa, D. Garcia, J. Guasch, R. A. Jimenez and J. Sola, Eur. Phys. J. C2, 
373-392 (1998), .hep-ph/9607485, 

[9] J. A. Coarasa, D. Garcia, J. Guasch, R. A. Jimenez and J. Sola, Phys. Lett. B425, 
329-336 (1998), |hipih7^11472| 

[10] J. Guasch and J. Sola, Phys. Lett. B416, 353-360 (1998), hep-ph/9707535j 

[11] A. Belyaev, D. Garcia, J. Guasch and J. Sola, Phys. Rev. D65, 031701 (2002), 
lhep-ph/0 105053} 

[12] A. Belyaev, D. Garcia, J. Guasch and J. Sola, JHEP 06, 059 (2002), |hep-ph/0203031[ 

[13] A. Belyaev , J. G uasch and J. Sola, Nucl. Phys. Proc. Suppl. 116, 296 (2003), 
|hep-ph/02 10253[ 

[14] M. Carena, H. E. Haber, H. E. Logan and S. Mrenna, Phys. Rev. D65, 055005 
(2002), hep-ph/0106116| [Erratum-ibid. D 65 (2002) 099902]. 

[15] A. M. Curiel, M. J. Herrero, D. Temes and J. F. De Troconiz, Phys. Rev. D65, 
075006 (2002), hep-ph/0106263 

[16] M. J. Herrero, S. Penaranda and D. Temes, Phys. Rev. D64, 115003 (2001), 
fhep-ph/01050^ 

[17] H. E. Haber, M. J. Herrero, H. E. Logan, S. Penaranda, S. Rigolin and D. Temes, 
Phys. Rev. D63, 055004 (2001), hep-ph/0007006 , 

[18] M. Guchait and S. Moretti, JHEP 01, 001 (2002), |hi^ph/0 110020 

[19] T. Plehn, Phys. Rev. D67, 014018 (2003), hep-ph/020612Tj 



10 



[20] T. Sjostrand et al., Comput. Phys. Commun. 135, 238-259 (2001), hep-ph/OOlOOlTl 

[21] A. Djouadi, J. Kalinowski and M. Spira, Comput. Phys. Commun. 108, 56-74 (1998), 
hep-ph/9704448, 

[22] S. R. Slabospitsky and L. Sonnenschein, Comput. Phys. Commun. 148, 87-102 
(2002), hep-ph/0201292, 

[23] E. Richter-Was, D. Froidevaux and L. Poggioh, ATLFAST 2.0 a fast simulation 
package for ATLAS, ATL-PHYS-98-131. 

[24] K. A. Assamagan, Y. Coadou and A. Deandrea, Eur. Phys. J. direct C4, 9 (2002), 
[hep-ph/02 03121( 

[25] K. A. Assamagan and Y. Coadou, Acta Phys. Polon. B33, 1347-1360 (2002). 

[26] H. L. Lai et al. (CTEQ Collaboration), Eur. Phys. J. C12, 375-392 (2000), 
lhep-ph/9903282, 

[27] S. Jadach, J. H. Kiihn and Z. Was, Comput. Phys. Commun. 64, 275-299 (1990). 

[28] M. Jezabek, Z. Was, S. Jadach and J. H. Kiihn, Comput. Phys. Commun. 70, 69-76 

(1992) . 

[29] S. Jadach, Z. Was, R. Decker and J. H. Kiihn, Comput. Phys. Commun. 76, 361-380 

(1993) . 

[30] D. P. Roy, Phys. Lett. B459, 607-614 (1999), |hSpph/9905542| 

[31] S. Moretti and D. P. Roy, Phys. Lett. B470, 209-214 (1999), |hepph/9909435 

[32] D. J. Miller, S. Moretti, D. P. Roy and W. J. Stirling, Phys. Rev. D61, 055011 
(2000), hep-ph/9906230 

[33] B. C. Allanach et al., Eur. Phys. J. C25, 113-123 (2002), hep-ph/0202233l [eConf 
C010630 (2001) P125]. 



11 



